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Abstract. We consider singular limits of the three-dimensional Ginzburg-Landau functional for a superconductor with 
thin-film geometry, in a constant external magnetic field. The superconducting domain has characteristic thickness on 
the scale e > 0, and we consider the simultaneous limit as the thickness e — ¥ and the Ginzburg-Landau parameter 
re — > oo. We assume that the applied field is strong (on the order of e _1 in magnitude) in its components tangential to 
the film domain, and of order log re in its dependence on re. We prove that the Ginzburg-Landau energy T-converges 
to an energy associated with a two-obstacle problem, posed on the planar domain which supports the thin film. The 
same limit is obtained regardless of the relationship between e and re in the limit. Two illustrative examples are 
presented, each of which demonstrating how the curvature of the film can induce the presence of both (positively 
oriented) vortices and (negatively oriented) antivortices coexisting in a global minimizer of the energy. 
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1. Introduction 

In this paper we continue the study of thin-film superconductors begun in our previous paper jA BGSlOj . 
The superconducting sample occupies a domain fi e CK 3 , 

n e = {(x',x 3 ) € M 3 : i'ew, ef(x') < x 3 < sg{x')}, 

where u C M 2 is a bounded regular domain in the plane, /, g : ui — > K are smooth functions on u> with 
f(x') < g(x') for all x' € to, and e > 0. We denote by 

a(x')=g(x')-f(x% 

the thickness of the film for given x' G ui. 

We study minimizers and Gamma-limits of the full three-dimensional Ginzburg-Landau model, for the 
superconductor f2 E subjected to a spatially constant external magnetic field, h ex G R 3 . The state of the 
superconductor is determined via a complex order parameter u : fi e C and the magnetic vector potential 
A : ]R 3 — > R 3 , which determines the magnetic field h = VxA. The energy of the configuration (u, A) is 
given by: 

I £ , K (u,A):=i f (|V^u| 2 + ^(l-|u| 2 ) 2 ) dx+i / |h-h c *| 2 dx, 
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with Ginzburg-Landau parameter k > and thickness parameter e > 0, as above. 

With the appearance of three parameters in this problem, e, k, and h ex , it is important to identify limiting 
regimes which are mathematically interesting and physically relevant. As in ABGS10J we rescale the domain 
by e in the x 3 direction in order to recognize the correct scaling for h ex in terms of the thickness parameter. 
We introduce rescaled quantities as follows: 

x = (x',x 3 ) = (x 1 ,x 2 ,x 3 ) = (xi,x 2 , £ Qi 

A{x) = (A 1 ,A 2 ,eA 3 )(x), 
u{x) = u(x). 

As a result, the order parameter u is defined in a fixed (e-independent) domain 

Q ■= O x = {(a:', a;) : f(x') < x 3 < g(x'), x' £ w}. 

We denote by h — V x A, and remark that the physical magnetic field h(x) = (e _1 /i'(a;), h 3 (x)) in the new 
coordinates. The energy transforms as follows: 

I e , K (u, A) = sI EiK (u,A), 

where 

l, K (u,A) =: J Q|(V - zA)u\ 2 + i _ ^ A 3 ) U | 2 + ^(1 - |u| 2 ) 2 ) d-T 

+^ 8 (i ft 3-ftri a +^K-fe' eiB r) 

and where V = ^2), and the rescaled effective external field takes the form 

h ex = (hl x ,h e 2 x , h ex ) = (eh?, ehf, h| x ) . 

In our previous work [ABGS10 , we considered the case where the Ginzburg-Landau parameter k was 
fixed, in the limit e — > 0. By the above transformation, we noted that taking the parallel component of 
the applied field h' ex — 0(e^ 1 ) gave a critical scaling for the applied field strength. We determined the T- 
limit of the energy in the case of critical, subcritical, and supercritical fields. The most interesting case was 
(unsurprisingly) the critical case. With critical parallel magnetic field, we observed an interaction between 
the parallel field component and the geometry of the domain, and the r-limiting functional was a two- 
dimensional Ginzburg— Landau-type energy, 

G K (v; F) = J^ a(x') {i|(V - *3>| 2 + ^(M 2 - I) 2 + M 2 } dx< , (1.1) 

with B' : M 2 — > M 2 satisfying V' x B' = F, an effective magnetic field F — F e 3 acting perpendicularly to 
the limiting plane of the film, 

F=h e 3 x - (h?, hf) ■ v ( ^mi±3^i ^ 

In this T-limit, the magnetic field is predetermined by the strength and direction of the original applied 
field h ex £ M 3 and by the geometry of the domain i\, and is part of the variational problem for the 
thin film limit. In this way, the functional obtained is of the same type as that studied by Chapman, Du, & 
Gunzburger CDG96 for thin film superconductors with constant k and bounded h ox . One of the attractions 
of the critical scaling case is that the effective magnetic field is non- constant when the film's vertical center 
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h(f( x ') + 9( x ')) i s curved, leading to some more interesting configurations of vortices (and antivortices) 
appearing in minimizers of G K near the lower critical field. 

In this paper, we consider the simultaneous limit of the energy as both e — > and k — > oo. We choose 
an exterior applied field which is critical with respect to the thickness parameter e, and on the scale of the 
first critical field in k, 

= (h'^H^k 

in other words, in the rescaled functional we take 

h ex — H log k, (1.2) 

where H — (H' , H 3 ) = (Hi,H 2 , H 3 ) € R 3 is a fixed constant vector (independent of e, k.) The choice of an 
applied field of order log k is natural for the Ginzburg-Landau model, both in two dimensions (see |SS07| and 
the references contained therein) and in three dimensions (see ABM06 , BJOS11J) as it is the critical scale 
of the magnetic field strength at which vortices become energetically favorable in the sample. Our results 



I e>K (u,A) := n _ I BjK (u, A) 



confirm this scaling in the thin film setting as well. For applied fields of the form (1.2 1, it is expected that 
the energy of minimizers of I e>K will be on the order of [logft] 2 . We are thus led to introduce the following 
normalization, and study the family of functionals 

1 

(log k) 

and configurations (u, A) with bounded values of I S , K - 

To present our results, we first introduce appropriate function spaces for the configurations (u, A). For u 
this is very simple, u 6 H (f2i, C). For A, we note that far from O, we expect the field h — curl A to relax 
to the rescaled field h ex . We first choose a fixed A with V x A = H = (Hi, H2,H 3 ). A convenient choice is: 

A := (h 2 x 3 - \h 3 x 2 , \h 3 Xi - Hix 3 , . (1.3) 

Note that this choice fixes a gauge for A. Then, we take our A in the following affine space, 

A e A := [A e Hl oc (R 3 ;R 3 ) : A - Alogrc G H\ iv (M 3 ;]R 3 )} , (1.4) 

where H^ iv (M 3 ; 1R 3 ) is the closure of the space of smooth, compactly supported divcrgcncc-frcc vector fields 
F e C^°(K 3 ;K 3 ) in the Dirichlet norm, (see }GP99j .) 

2 



\DF\ 2 dx 



With the onset of vorticity, the limiting behavior of the functional I E K must be described in terms of 
the limiting currents and vorticity measure (Jacobian) rather than the order parameter u, as the number of 
vortices will become unbounded in the limit. The current (or momentum density) is: 

3 — j( u ) = (iu,Wu), where (a, b) = - [ab + ab], 

and the vorticity, or weak Jacobian, J = |V X j. For u S _ff 1 (il;C), j £ _L 2 (0;M 3 ), and so J is defined in 
the sense of distributions, although in our context it will in fact be measure-valued (see Jerrard & Soner 
[JS02 .) It will often be convenient to represent j and J as differential forms, 

j = j! dx 1 + J2 dx 2 + ]3 dx 3 e A 1 ^ 3 ), 
J = dj = Ji dx 2 A dx 3 + J 2 dx 3 A dx 1 + J 3 dx 1 A dx 2 e A 2 (K 3 ), 
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as the natural mapping between forms and vector fields is an isometry in Euclidean space R 3 . The F-limit 
will impose further structure on the limiting momenta and Jacobians, so we define the following domain: 

Z := jj G L 2 (n-R 3 ) : j = (j'(x'),0), J := x j G 2tt(0;IR 3 ) j , (1.5) 

where 9Jt(f2; K 3 ) is the space of vector- valued Radon measures on fi. Note that for j e Z, the corresponding 
Jacobian takes the form J = (0,0, J 3 (a/)). We define the functional 

i-\\a(x')V x jllmtiv, + - / a(x') \f - B'\ 2 , if j G Z, 

iooU;F) = \ 2 zJu (1.6) 

[oo, otherwise, 

where F = V' x B', and B' : R 2 — > R 2 . This is our mean field model, or vortex density functional (see 
|CRS96| .l 

The main result is that 7 £jK Gamma-convergences to Ioq. We prove this in the usual two steps: first, 
bounded sequences are compact and the limit is lower semicontinuous in the energies: 

Theorem 1.1. For any pair of sequences e n — > and n n — > oo, assume {(u„, A n )} nl z^ C i? 1 (17;C) x A 
satisfy the uniform bound, 

sup I eniKn (u n) A n ) < +oo, 

ngN 

and define j n — j(u n ) = (u n ,du n ) and J n = hdj n . Then there exists a subsequence (which we continue to 
denote {e„, n n } ) and G Z , with = ~V x j*, such that: 

(1) along the subsequence, 



| 2 ^1, mi 4 (fi)„ (1.7) 
A^O, in H\ iv {R 3 ;R 3 ), (1.8) 

j», mL 2 (n ; R 3 ), (1.9) 



log ft 

3 n 



l0gK 



| Ural logK 

J + , m t/ie weak--k topology on [C°' 7 (f2)]*, (1-10) 



for all < 7 < 1. 
(2) Furthermore, 



1 /" a 3 (a;') ,2 
liminf 7 e ljK „(w„, A„) > Ioc{j*;F«) + - / — — — Lff' dx', 
n-s-oo 2 J w 12 



where is defined by 



F*{x') =H S - (H l7 H 2 ) ■ V f /Ml+ffM ') ^ (Ln) 



and Ioc(j\ F) is defined as in (1.6). 

It is important to note the (somewhat surprising) fact that the same limit is obtained regardless of how the 
two parameters e„ — > and k„ — > 00. In fact, this is particular to the case where h ex = 0(log k). For stronger 
applied fields h ex 3> log k the character of minimizers will depend on the relationships between k, e, and 
h ex . In particular, in Remark pO| we note that when e is relatively large compared to k the superconductor 
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will not behave as a thin film at all, and may exhibit a longitudinal vortex lattice, aligned along a horizontal 
direction. 

The proof of Theorem |1.1| is the content of section 2. The second part of the Gamma convergence result 
is the construction of recovery sequences: 

Theorem 1.2 Let j G Z and consider any sequences e n ,K n such that e„ — > and n n — > oo. Then there 
exists a sequence {(u n ,A n )} C i? 1 (r2;C) x A. satisfying 



J„ 



log K n 

Jn t 1 



j in L p (fl), for all p < 2, 



J := -V x j weakly in M(il;R ), and strongly in (CJ(fi))', < 7 < 1, 



log K n 2 
with j n := (iu n ,du n ) and J n :— \dj n . Moreover, 

liminf I Sn>Kn {u„, A n ) < i 0o (j*;i ? *)- 



We will prove Theorem |1.2| in section 3. The construction is essentially two-dimensional, but since the 
effective applied magnetic field — V x B^ is non-constant the procedure is somewhat different from the 
standard approaches with a constant applied field. 

An immediate consequence of the Gamma convergence of I e . K to 1^ is the convergence of minimizers. To 
better understand the properties of minimizers of the mean-field limit 1^ we use convex duality to obtain 
an equivalent formulation of the problem as a variational inequality. 

Proposition 1.3. Assume j £ Z is a global minimizer of 1^, and : M 2 — > K 2 with Vxfi, = F*. Then 
j = B. t — ^V^C where ( G Hq(lj) solves the minimization problem: 

1 



min £oo(0, with £00 (0 := / 

= H <U<") Jul 



iei<»(=)/2 



a{x) 



(1.12) 



The equation for C is a two-obstacle problem for Poisson's equation, and leads to solutions with free 
boundaries on the coincidence sets, where |CI — ^a* ^-he coincidence sets form the support of the measure 
J, and indicate the regions of nonzero vorticity of minimizers in the simultaneous K — > 00, e — > limit. 
Obstacle problems of a similar type were obtained in the limit of the two-dimensional Ginzburg-Landau 
functional with constant vertical applied field by Sandier & Serfaty ISSOOj . and in a non- homogeneous 
setting involving pinning of vortices in [ASSOlj . In these papers, there is a single obstacle, as the analogue 
of the solution £ is constrained on only one side, and in the 2D setting the equation is of Hclmholtz (rather 
than Poisson) type. 

This more concrete characterization of the limiting problem gives us a better idea of what minimizers 
look like for applied fields on the order of the first critical field. In section 4 we revisit two examples (for 
superconducting films which approach a disk in the e — > limit,) which we introduced in our first paper. In 
the first example, vortices accumulate in two symmetrically placed subdomains in the disk, one containing 
positively oriented vortices, and the other antivortices (with negative winding.) In particular, this implies 
the rather surprising conclusion that vortices and antivortices can coexist in global minimizers of the three- 
dimensional Ginzburg-Landau model with a constant applied field. 

The second example illustrates the phenomenon of concentration on curves and annular subdomains, 
which occurs even in a simply connected domain in the thin- film limit. This example is quite appealing in 
that the free boundary problem obtained is radially symmetric, and an explicit solution may be calculated, 
and the changing geometry and topology of the regions of vorticity are explicitly shown. Previous examples 
of concentration on curves or annuli were found for minimizers in annular domains (see |AAB05| . [AB06 , 
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[AB05J, [ KaclOj . |Rou| .) Determining the asymptotic distribution of vortices on curves in the domain requires 
more delicate estimates very near to the lower critical field; this is done in ABM11 . 

Finally, a different type of thin film problem has recently been studied by Contreras & Sternberg CS10 
and Contreras [Conllj . In their setting, the superconductor is thin shell, built from depositing an e-thick 
coating on a fixed two-dimensional surface in IR 3 . The limiting problem in this case is a Ginzburg-Landau 
model on an embedded 2-manifold, and they obtain remarkable results connecting the lower critical field 
and the appearance of vortices to the geometry of the limiting surface. 

Several results on Gamma convergence or the convergence of local minimizcrs have also been proven for 
three-dimensional models of superconductivity (or Bose-Einstein condensation) without assuming thin-film 
geometry. A recent paper by Baldo, Jerrard, Orlandi, & Soner [BJOS11 proves Gamma convergence under- 
very general hypotheses. Results on stable vortex solutions for these models in various specific geometries 
may be found in [MSZ04] . |.IMS04j . [Jer07] . [ABMOBj . |ABJV108j . 



2. Compactness and lower bound 



In this section we prove two parts of the Gamma-convergence result, the compactness of energy-bounded 
sequences and the lower bound inequality. As always when dealing with the Ginzburg-Landau functionals, 
the gauge invariance of the functionals is an issue. The choice of the space A fixes a gauge (see (1.3) and 
(1.4) above.) We will require the following essential lemma: 

Lemma 2.1 (Lemma 3.1 in |GP99j K Let g £ L 2 (M 3 ;R 3 ) such that divg = in X>'(IR 3 ). Then there is a 
unique B E H 1 (M 3 ;M 3 ) such that V x B = g and div B = 0. 



As a consequence, it follows that 



\B\ 



IV x B\ 2 dx 



is equivalent to the usual (Dirichlet) norm on the space H\ iv (1 



5 ). 



Proof of Theorem 11.11 

Let K := sup n g[ij Ie nt K„ (Wjj, j4«) < oo. From the energy bound we immediately deduce that, 



(Kl-i) 4 < / (\u n \ 2 -iy< 



, K(\nK n ) 2 



hi 



log K v 



H' -> in L\ 



K n 

2 / TTD 3 . TTT) 2 \ 



0, 



(2.1) 
(2.2) 



In particular, |u„| — > 1 in L 4 (fl). The vertical components of the magnetic field are bounded via the energy 
bound, and thus along a subsequence (which we continue to denote u n ,A n ), we may conclude the weak 
convergence, 

h 3r 



H] 



log K n 

H := (0,0, 1) in L 2 (K 3 ; 



in L Z (I 



(2.3) 



s ), and each div h n = (in the sense of distributions,) 



As the vectors f , hn 

L log K u 

we may conclude that the limit is also divergence-free, div if = 0. As a consequence of (2.3) and Lemma 

! ; M 3 ) with V x A = H = (0, 0, t) and 

A 



we also conclude that there exists A g H^ iv 



2.1 
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weakly in H j iv (K 3 ; R 3 ) , and in the norm topology on LP(Q), 1 < p < 6. Since H = (0,0, £) € L 2 (K 3 ;K 3 ) 
with div H — 0, we conclude that = (in the sense of distributions,) and thus I — 0, and also H = = A 
(by Lemma 2.1 ) In particular, (2.4 1 implies 



log K„ 



A — > weakly in H^, 



p3. ttc>3' 



, and in the norm on L P (S1), 1 < p < 6. 



(2.5) 



To obtain the lower bound we adapt the argument of [SS04 . Expanding the quadratic term in the energy 
bound, we obtain: 

2K>2I EntKn (u n ,A n ) 

^ (log«„)- 2 f (|W| 2 - 2A' n ■ (iu n ,Vu n ) + M 2 K| 2 ) dx 



> 



n 



1 




2 


A' 


2 


log K n 




log K n 


1 


V'u n 


2 


4' 


2 


log K n 




log K n 



\u n \ 2 -l)- 



A{ 



log K T . 



By (2.5 1, the last term is bounded, and 

A' 



<> 



log K r 



1 < 



log K r 



I |2 ill / ^ lo g K ™ 
\Un\ - l\\ L 2 < C . 



K I, 



A' 



by the energy bound and the L p boundedness of t ^ . Thus we have 



< c, 



1 




In 


log K n 



(2.6) 



with constant C depending on the energy bound K. By a similar calculation, we may also obtain the estimate 

2 

dx < C, 



[1 


d 3 Un 


In £ n 


log K n 



and so we have strong convergence in the a^-direction 



inZ/(Q;C) 



log K n 

We now turn to the currents, j n := (iu n , Vu n ). Following |JS02] we normalize the currents as follows 

~ jn _ (iu n ,Vu n ) 

\u n \ \ogK n \u n \ log 

We observe that each component of j n = (ii, n , J2,n, J3,n) is (for fixed n) pointwise (a.e.) bounded, 



(2.7) 



,~ I , \d k U n \ , 

Ufe,n < 1 , « = 1,2,3, 

log K„ 



(2.8) 



so jn is well defined almost everywhere in f2. Moreover, from (2.6) and (2.7) it follows that there exists 
3 = (i'jO) € -L 2 (£!;R 3 ) such that (along a subsequence) 



in i 2 (S7). Writing 



in ^ j"i i3,n -> 



log K 
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we recall that (\u n \ — 1) — > in L A {Vt) (see (2.1 ),) and thus obtain that 



J3,n 



-> 



L 4 / 3 (n). 



(2.9) 



log K n log K r , 

We will require this fact later on in determining the limit functional. 

We continue as in the proof of Theorem 2 of [SS04] . with some simplification due to the thin film limit, 
and our modification of the currents. Let ei,e2,e3 be the standard basis in M 3 , and define vector fields 



X k = h k e k , k = 1, 2,3, with h k £ C (ft) and |/i fc (x)| < 1 for all x e Q, k = 1, 2,3. By ( pL6[ ), we have 

|X fe -V'u„| |/lfe9feM n | 



log Kn 



log K„ 



^x fc , fc = l,2,3, 



(2.10) 



weakly in L 2 (f2). Passing to the weak limit in the bound (2.8), we conclude that 



\X k -j\ = \h k j k \<<f> Xk , fc = l,2,3, 
pointwise a.e. in f2. By (2.7), <f>x 3 = 0. We also define the defect measures vx k corresponding to the weak 



convergence in (2.10) 



X k ■ V'ttr, 



log K n 



+ i^x fe in the sense of measures, 



(2.11) 



for k — 1, 2,3. Because of the strong convergence in (2.7), it follows that ^3 = 



For the Jacobians J n , we apply Theorem 1 in [SS04j : by the energy bound, 



H\u n \ 2 -l) 2 ) <C[ln Kn ] 2 



with constant C independent of n (using the estimates (2.6), (2.7), and (2.1)), we may conclude that 



log Kn 



± J* , 



in the weak* topology on [C ' Q (f2)]*, for < a < 1. Moreover, the limiting Jacobian is a Radon measure- 
valued two-form. Furthermore, the same theorem relates the limiting Jacobian to the defect measure vx k 



via a product formula (see (2.12) below.) We prove the following properties of the limiting Jacobian and 
currents: 

Lemma 2.2. The limiting Jacobian J* = \dj* has the form J* — J| dx 1 A dx 2 with J3 = J^{x'), and the 
limiting current j* G 2, . 



Proof. We make use of the product formula from |SS04| in the case where N K = 0(ln«;), which we review 
here. Let £ be a bounded smooth domain in K 3 , and v K £ H 1 (E,C) satisfying 

E K (v K ;E) := jf (^|V« K | 2 + ^(\v K \ 2 l) 2 ) < C[\n K } 2 , 

for constant C independent of k. Let X, Y be continuous, compactly supported vector fields in E, and ^x, 
vy the defect measures (defined as in (2.11)) for v K as k — > 00. Then, the normalized Jacobians — v J in 
(Cq(Q,))' for all 7 > 0, and the defect measures are related to the limiting Jacobian via: 

2 



Wx\(E)\v Y \(E)> 



J {XX) 



(2.12) 



Here we denote by H(-E) the total variation of the measure v over the set E. 
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We note as above that for any E C 12, E Kn (u n ;E) < [In K n ] 2 / e „ lK „ (u n , A n )) < C[ln/c„] 2 . Let E be any 
open ball contained in fl, and Xk — hk&ki with hk G Cq{E) and |/ife| < 1, fc = 1, 2, 3. Applying the product 
formula we then obtain, 



0=\v Xl \HE)\v X3 \HE)> 



hi h 3 J*(ei,e 3 ) 



Taking the supremum over all such hi, h 3 , we conclude that, as a Radon measure, J*(ei, 63) = in the ball 
L. By an analogous computation with X 2 = h% e 2 and X 3 (as above), we also have J*(e 2 , 63) = in the ball 
E. This holds for any ball E C fi, and thus these measures vanish identically in f2, and thus the Jacobian 
has the form J* = J3 dx 1 A dx 2 . 

Furthermore, since J n — \dj n for each n, it follows that dJ n = (in the sense of distributions.) Normal- 
izing by logK„ and passing to the limit, we retain dJ* = 0, and hence d 3 J^ = in D'(fi), so Jj* = J£(x'). 
This also implies that (in the sense of T>'(VL),) 

= J* = 9 2 i*3 - 9ai*2 = -<9 3 j*2, 
= J 2 * = 93j*i - 9ii» 3 = 9 3 j*i- 

Thus, the limiting current must have the form j t = (jl(x'), 0) and J* — gV x J* S 9Jt(f2), and hence j* E Z. 



It remains to verify the lower bound inequality. From the definition of the defect measures and the 
product formula from Theorem 1 of |SS04| . 



lim inf 



/ 




In 


log K n 



> V lim inf / 



X k ■ V'it„ 



fc=l,2 

> \v Xl \(Sl) + \v X2 \(Sl) 



> 2 
= 2 



log Kn 

/ J*(X!,x 2 ) + / (AVj*) 2 + (* 2 -j*) 2 

Jn Jn 



hi h 2 J*(ei,e 2 ) 



+ / ei | 2 +^|i*-e 2 1 2 ) 



(2.13) 



The above estimate is valid for any /i^. € Co(fi) with |/ij.(a;)| < 1, k = 1,2. We choose these functions to 
obtain an estimate in terms of the total variation of the measure J 3 = J*(e l5 e 2 ). By the Hahn decomposition, 
we may write J3 = /i + — /i_ for mutually singular, nonnegative finite measures supported on the 

disjoint sets E + ,E_ e f2, respectively. Take sequences hi^,h 2 .i £ C (Q) with < 1, fc = 1,2, such that 

hi ti -t 1, h 2 ,i -> xb+ - X-B- 

pointwise a.e. in 17. Passing to the limit i — > 00 on the right hand side of ( 2.13[ ) (using the Lebesgue dominated 
convergence theorem,) we conclude that 



lim inf 

n— > 00 



/ 


V'u n 


In 


log K n 



>2\j;m)+ [ 

Jn 



(2.14) 



Finally, we derive the form of the lower bound for the full energy. First, from the strong L 4 convergence 



(2.1 1 of \u n \, the weak La convergence (2.9) of the normalized currents, the strong L p (1 < p < 6) convergence 
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of the vector potentials (2.5), and the lower bound (2.14), we may conclude that 

1 

n— f oo 



liminf I Sn ,«„(w„,A0 ^ Umirif 

n— >-oo 

> IU*IU 



2[log« 1~ 



(I W| 2 - 2A' n ■ f n + K| 2 + (|u„| 2 - 1) K| 2 ) dx 



raj "'a 

|j;| 2 -2A.j; + |1| 2 



Since both J 3 = J3(ie') and — j'*{x'), we may integrate out the variable x 3 , to reduce to a two- 
dimensional total variation, weighted by the film thickness function a(x'), 

H-^HsuKn) = Il a (^ )^3||an(aj)- 



The limiting vector potential A (defined in ( 1.3 )) is x 3 -dependent, but this dependence may be averaged out 
(to produce the desired effective field F*). Indeed, we decompose A' as follows: 



#32:2 , ^#3^1 ) + (H 2 x 3 , -Hxxz) =: A± + (H 2 x 3 , -H x x 3 ) 



A' 

Expanding the energy and integrating out x 3l we have: 

\j:-A'\ 2 dx= f \ J :-A ± ~(H 2 ,-H 1 )x 3 \ 2 dx 



a{x')\ji - A ± \ 2 - 2 / (jl - A±) ■ (H 2 , -H x ) / x 3 dx 3 dx' 



g( x ') 



/(*') 



\H 2 , —Hi) 



x\ dx 3 dx' 



a(x'M-A ± \ 2 -2 I atf)£l&(ji-A ± ).(H 3 ,-H 1 )d a > 



a(x') 



>\ (r+/g+g 2 ) 
3 



(H 2 ,-Hi)\ dx' 
2 



^l\H'\ 2 dx'. 



We conclude that 



= J otf)\jl-A ± - {^{Hz-Hx 
lnninf7 En , K „K, A n ) > \\a{x')J 3 \\ m(u)) + a(x') (\£ - B'f + ^1\H'\ 2 ^ dx' 



with 



B' := A ± 



f+g 
2 



(H 2 , —Hi) 



Since V' x B 1 ^ = with as given in ( |1.11[ ), this concludes the proof. 



(2.15) 



(2.16) 



Remark 2.3. The fact that the same T-limit is obtained regardless of how the two parameters e n — > and 
K n —¥ 00 is somewhat remarkable. In fact, this is particular to the case where h ex = 0(\ogn). To see this, 
take a rectangular solid domain Q £ = (0,a) x (0, b) x (0,e) with horizontal applied field h ex = ( , 0, 0). 
The two-dimensional minimizer of the Ginzburg-Landau energy (w K (x 2 , x 3 ), A K (x 2 , x 3 )) may be used as a 
test function, with energy (see Theorem 8.1 of |SS07j ) 

A K ) * \n e \ f log " * hm log J?l- 
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A "thin film" solution, that is, one with u = u(xx,X2) and A — A has energy given by (1.1 1, with effective 
field F = 0, that is 



:Gju) = : / uu')\ i|V U | 2 + ^(\u\ 2 - I) 2 + ^l\ Hl ( K )\ 2 \u\ 2 } ,/.,'. 



For re large, the minimizer is essentially u = 1, and so a thin film configuration has energy of the order of 
eH^n) 2 . Thus, if 



2 

re z e 



Ht(n) \og--^— < £#i(re) 2 , 

we expect that the minimizers do not have thin film form. To take a concrete example, if the magnetic field 
Hi{k) = re, then the horizontal vortex lattice is preferable for e 3> - with re large. 



3. The recovery sequence 

This section is devoted to proving Theorem |1.2| namely the existence of a recovery sequence and the T — 
lim sup inequality. In both this section and the following one we require the following Hodge decomposition 
with respect to the weighted inner product, 



(v, w) = j a{x)v-wdx' 
on L 2 (u;;R 2 ). We define the following subspaces: 



v = {vc, cefrVjR)}, (3.1) 

W = {W £ C\uj;R 2 ), V- 1 • W = 0, V • (aW) = 0, W ■ v = on du>) . 

Lemma 3.1. Any Z G L 2 (lo;M. 2 ) admits a unique orthogonal decomposition Z = U + V + W with U £14, 
V G V, W G W, with respect to the inner product (■, ■). The space W is finite dimensional: in case uj is 
simply connected, W = {0}, and m case uj — ujq \ U^LjWj /ias m holes, dim(W) = m. 

Proof. First, assume Z G C°°(a;;IR 2 ). We define ip and £ as the solutions to the boundary- value problems, 
_V • ( -^Vip) = curlZ in w, [ V ' « x )^0 = div in w ' 

^ = o on ftj, I ^ = z - v 011 9w > 

Then, it is easy to verify that W :— Z + ^ ±r ^ — VC satisfies curlW = = div [aW] in us, and W ■ v = 
on duj. Moreover, by integration by parts we see that W _L i-V ip -L V£ in the inner product (■,■). 

To identify the space W, we apply Lemma 1.1 of BBH94 and note that any W G W may be written 
as W = ^V -1 "^ with £ constant on each component of du>, and V • ^V£ = in u. If a; is simply connected, 
the maximum principle ensures that £ is constant in uj, and W = {0} is trivial. In case uj — uj \ U™ x uji 
is multiply connected, we follow the treatment of AB06 . For each fixed i = 1, ...,m we define functions 
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£i G Hq(lu) which solve 



V • -V& = 0, in ui, 
a 

€i\duj = <kj, j = l,..., m 
' ^ dx = Si j j = 1, . . . ,m, 



2"" Jdu, a dv 



(3.2) 



where c^- are constants (determined by the solutions,) and 5i.j is Kronecker's delta. The existence of such 
£i may be obtained by minimizing 

m) = l I -|Ve| 2 dx + 27reL 1 
2 Ju a 

over the class of £ G i?Q (wq) with £| w constant. (See section 1.1 of [BBH94I .) It is easy to show that 



Thus, W = £ G W is parametrized by the m constants $i, i = 1, . . . , m, and W is m-dimensional. By 
elliptic regularity we also have W C C 1 (w;IR 2 ). 

For general Z G L 2 (w;R 2 ), the general result is obtained by density. 



We are now ready to complete the proof of the Gamma convergence result. 



Proof of Theorem 1.2 



Let j G Z be given, as well as a sequence n n 



oo. We choose vector potentials 

^4 e:E , and construct a sequence of order parameters u n of the form u n (x) = v n (x') to satisfy the demands 



of the theorem. As noted in [ABGS10 , for configurations of this form, the three-dimensional energy reduces, 

I Sn , Kn (u n ,A n ) = G K (v n ;F»), 



where G K is defined in (1.1) and is as in (2.16). Since all which follows will be two-dimensional, we drop 
the primes in our notation, and write -B* G K 2 , V = (d Xll d X2 ), and v{x),a{x) for x — (x\,X2) G K 2 . 
We next apply the Hodge decomposition above to our given j G Z, and write 



j = U + V + W 



1. 



-V x ip + vc + w, 



with -0 G Hq(uj), ( G H 1 ^) and W G W, a mutually orthogonal splitting in the inner product (•,•). Since 
V, W are irrotational, they do not contribute to the weak Jacobian J = x j, and carry no vorticity. As 
in |JS02| . we may associate to V, W an S^-valued map w K . The singular part of the Jacobian is contained 
in U; for this part we construct a family u K with point vortices via an appropriate Green's function. We 
adapt the arguments of [SSOOj to deal with the inhomogeneity of the functional /. Putting these two parts 
together, the desired recovery sequence will have the form v n — w Kn u Kn . 
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Step 1: The components V + W e V © W 
From the proof of Lemma 

E£i$<6(a0> for 6 a s in © 
denote the integer part, and set 



From the proof of Lemma 3.1 we may write V = V£, ( € H 1 ^) and W = ^V -1 "^ with = 
for fa as in (3.2) with real constants. Let M i: „ = [$jln/e n ], i = 1, . . . , m, where brackets 



S n := V M itn d, W n = — W^E n . 
U 

We note that 

\\W n -W\nK n \\ c , <C, (3.3) 

for constant C depending on W (but independent of n.) 
Since 

curl W n = V Mi.nV- 1 • iv x & = 0, and / VK„ • r ds = V M hn I - ^ds = 2 7 rM i , n , 

an integer multiple of 2ir for each j = l,...,m, it follows that W n is locally a gradient, W n = V?7„ for 
77 n possibly multiple valued, but for which e lv ™ is smooth and single- valued in lj. We may then define the 
complex order parameter 

w n = expi(r) n + (lnn n ). 

By construction, 

iM J w - Vw "> 4 nff (3.4) 

log K n log Kn 

in C 1 (w). Since |w n | = 1, we may easily calculate the contribution to the energy using the orthogonality: 

a(x)\\7w n \ 2 dx — — / a(x)\\7rj n + V( In n n \ 2 dx 



2./,., v " 2 



' f a(x)\W n \ 2 + ^^- I a(x)\WC\ 2 dx 



2 



< / a(*){|Wf + |l'| 2 }da ; + 0(l) ) (3.5) 



using (3.3) in the last line. This completes Step 1. 

The treatment of the component U = — ^V^V- 1 £ W will require several steps. First, we restrict to 
ip G Cg°(u>); the result for general ip £ Hq(uj) will follow from a diagonal argument. Denote byKdw the 
support of ip. 

Step 2: Approximating the measure fi := cur\U = — V x ;jV _L ?/> by Dirac masses (representing vortices.) 
Let N n S N be any sequence of whole numbers with 
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1 N n 



i=l 



Applying Lemma 7.5 of [JS02 ], there exist families of points {p™}i=i,...,jv ra m the set X = supp-0 and 
associated integers cr™ G { — 1, 1} with the following properties: 

\Pi ~Pj\> cqN~ 1/2 for i ^ j, for constant c = c o (-0); (3.6) 
limi?(a) = where R(a) = limsup V ^ (3.7) 

\pf -P™ I 

Mn := aT a i S p" ~^ V, (3-8) 

n i—1 

1^1 = ArE<W^H> (3.9) 

n i—l 

where the convergence in (3. 8), (3. 9) is weakly in the sense of measures, and strongly in [Cq' 7 ]' for all 
< 7 < 1. By |ju| we mean the total variation of the measure fi = curl U. 

As in [SSQO] we modify the measures /i„ by regularizing the Dirac mass. Let /i™ := KnH 1 ldB(pv,i/ Kn )i 
the element of arclength on 5" :— dB(p™, l/n n ), normalized with mass 2ir. We define the measures 

withp" € K, erf G {0, 1} as above. Since each /i™ — ► 5 p ™ strongly in [Co' 7 (w)]' for all < 7 < 1, and weakly 
in 2t(w), we may conclude that (3.8 1,(3.9 1 hold as well for u n , 

v n — > fx, \v n \ — > strongly in [Co' 7 (w)]' and weakly in 9Jt(w). (3.10) 

By Fubini's theorem we also note that the product measures also converge, 

v n ®v n — (3-11) 

strongly in [Cq' 7 (w x uj)]' and weakly in 9Jt(c<j x uj). 

Step 3: Recovering U from /1 = curl U. 

We introduce the Dirichlet Green's function, G a (x,y) in uj, which solves 

-Va, • —r -rV x G a {x,y) = 6 y (x), in u, 
a(x) 

G a (-,y)=0, on dui, 

for each fixed y G uj. By standard elliptic theory (recall a > is smooth in uj) we may conclude that G a (x, y) 
is smooth in uj x uj \ {y — x}, and 

G a (x, y ) = -^ln\x-y\ +1 (x,y), (3.12) 
where the regular part 7 has the property that for every compact set K <s ui, there exists C(K) < 00 with 

sup \j(x,y)\ < C(K). 

veK 

■ EE 

Given U G U, we then obtain the potential function ip G Hq(ui) from curlC/ = fi by solving 

I ip = on 9cj, 
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and we recover U = -\V^ip. Using the Green's function representation, we have 

ip(x) = j G a (x,y)dfj,(y). 

J u 

Since /i £ H~ l (ui), we may calculate the weighted norm of U in terms of the measure [i as follows: 

a(x) \U\ 2 dx = [ - IV^I 2 dx 



t/> ■ I -v-> ) dx 

ip(x) d/i(x) 

G a (x,y)dn{y) dfi(x). 



(3.13) 



Step 4: There exists a sequence ip n £ Hq(lo) for which — -V Vn — > U strongly in L p (uj) for all p < 2, and 

1, 



lim sup / ~\Vip n \ J dx< j a(x)d\(i\(x) + / a(x)\U\ 2 dx. 



(3.14) 



For each n, we define ipn{x) — j G a (x, y) dv n {y), and so ipn solves 



on 9a; . 



By (3.10) and elliptic regularity, we have ipn — > ^ in W 1,p (oj) for allp < 2, and thus — rV ^ n — > U in L p {u>) 
for all p < 2 as claimed. 

To estimate the energy we use the Green's representation. Since v n £ i? _1 (a;) for fixed n, by (3.13) we 
conclude that 

1 



\Vlpn\ dx 



G a (x,y) dv n (y) dv n {x). 



For any < a < 1, let A a = {(x, y) £ lu x lu : \x — y\ ^ a}. Fix \a G C°°(cj x w) with < % a < 1, and 

1, if a; £ A a , 



0, ifx^A 2ct . 



For any a £ (0, 1), G a (x, y)(l — Xa(x, y)) is smooth, and hence by the strong [Cq' 7 ]' convergence v n — > \jl we 
have: 



lim 

n— ^oo 



G a (x, y)(l - Xa(x, y))dv n {y) dv n {x) 



Xa (x, y))dfi(y)dn(x). (3.15) 



For the complementary integral, we use (3.12) to observe that 

a(x) , ] 



G a (x, y)Xa{x, y)dv n (y) dv n {x) 



K JA 



log 



2tt \x-y\ 



l{x,y) 



dv n (y) dv n (x) 



< 



K JA 2a 27r 

1 Nn 
J2 51 



a(x) 1 

log — — — dv n (y) dv n (x) + Ca 



\x - y\ 

°m log dfl7{y) d7{x) + Ca 

2vr \x-y\ 



(3.16) 
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To evaluate the remaining integral, we consider the contribution due to distinct points pf ^ in A2 a 
separately. We adapt an argument in Proposition 7.4 of [SS07 . Define the index set 

J n = {(i,j): \ P 7-P r l\<2a}. 



Let R n = \cqNu ' , where cq = Co(V>) is the constant in (3.6). We also define balls S™ = B(pf,R n ), 
i = 1, . . . , N n . By the choice of they are disjoint, as is the union 



(J (Bi x Bjj C A 3Q . 

We also observe that for any R < R n and (i, j) <G J n , since R < \\Pi — Pj\, we have 

\ < < ^ for aU x g B(p ?, y g fi(p n 

For (i, j) G J n we then have (recalling that = dB(pf, = supp/x™,) 



log | 1 

B"xB" F i/| 



da; dy > 



B"xB" IP, P 7 ' I 



> 



4 

J2_ 



|J>?-P?I 



lQ S |„n „n 

9,"xS™ \Pi Pj 

log I 1 



■dfi{x)dtf(v), 



(3.17) 



using (3.171 in the first and last lines. Summing over all pairs € and using the disjointncss of the 
union of the Bf x B™, we obtain: 



L y 



iV2 



S"xS' 



log I .da; 

B"xB" F — yi 



< c 



log, rdxdy =: 11(a). (3.18) 

\x-y\ 



As | log \x — y\\ is integrable, the remainder 1Z(a) — > as a — > 0, and so this term will not contribute to the 
limiting energy. 

Finally, we consider the contribution from the self-energy of the vortices p™. We parametrize the integrals 
over S?* = dB(p™, using complex notation, that is we write x,y £ dB(p™, —) as x = p" 
y = + ±~e iT , 0<6,T< 2ir. Then we have: 



1 AS 



1 

JV2 



a !^log v ^—dtf(y)dtf(x) = -^ 
2tt \x - y\ N£ 

1 

~ K 
1 



2, ^ a (p» + £) 



JO 

2ir 



2tt 



log K„ + log 



i(B- T ) _ 1 



d6 dr 



a\Pi + 



d9 + 0{N n 



a(x)d\tf\(x) + 0(N-- 
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Summing over all i = 1, . . . , N n , we arrive at 

iv2 



a(x) 



1 



log, r dM?(l/)dM?(a;) 



2tt 



a^d^lM + OliV- 1 ). 



(3.19) 



Passing to the limit n n — > oo, we thus obtain from (3.15l,(3.16),(3.18), and (3.19), that 



limsup / / G a {x,y)dv n (y)dv n {x) 

< / a{x)d\n\{x)+ / / G a (x,y){l-x a (x,y))dfi(y)d^{x)+Ca + Cn{a). 

J U) J UJ J u 

By hypothesis, the measure /i is bounded and absolutely continuous, and so we may apply dominated 



convergence to pass to the limit a -> and obtain the desired bound (3.14), as 



limsup I -\Vijj n \ 2 dx = limsup / I G a {x,y)dv n (y)dv n {x) 

< / <z(x) d\n\(x) + / / G a (x,y)dn{y)dn(x) 

J CO J UJ J UJ 

a(x)d\fi\(x) + / a(x) \U\ 2 dx, 



by (3.13) 



Step 5: Construction of a sequence u n £ Hq(lj;C). 

Let U n = -iV n ±V- L Vn- Then, V^Un = N n V ■ (^V?/>„) = locally in w \ uf"B(pf , i). Moreover, if C 
is a simple closed curve inw\ U 1 ^ rl B(p", we have 



U n • t ds £ 2-7rZ, 



by the normalization = 27r. Thus, we may write C/„ = V</>„ in oj \ U^B^pf, with <fi n which is 

multiple valued, but for which V</>„ and e"^ n are single- valued inw\ L)f n B(pf, -^-). 
To remove the singularity at each vortex core we define, 







if i*-pri<^ 



p?{x) := <j 2K„|a! -p?| - 1 if ^ < \x < £ 
1 if|z-p?|>i 



and ,o„ := J^J /?". A simple computation shows that 



i=l 



"i-i^ivpri 2 + ^((p- i ) 2 -i) 2 )^».- o,. 



with constant Cq independent of n. 



18 S. Alama & L. Bronsard & B. Galvao-Sousa 



Now define u n — p n e , with p n , (f> n as in the preceding paragraphs. We then have: 

{^|Vu„| 2 + ^(K| 2 - 1) 2 | dx = J a(x) {~p 2 |V0„| 2 + i|V P „| 2 + ~{p 2 n l) 2 | 

/V 2 r i 

2 A, a(ac) 



From (3.14) we then conclude that 

limsup— f a (^)(i|V Un | 2 + 4(kn| 2 -l) 2 )^< \ [ a(x)d\n\(x) + l [ a{x)\U\ 2 dx. (3.20) 
Since (p 2 — 1) — > in L q for all g < oo, we also conclude that 

= --V X Vn + ^ ~ Pn) V^^n — » f7 in for all p < 2. (3.21) 

iV n a a 

Step 6: Putting it all together. 

This follows as in [JS02]; we provide details for completeness. Write j £ Z as j = U + W with U £1A 
and W = V + W, V £ V, W £ W. Let w n be as defined in Step 1 and u n as constructed in Step 5, and 
define v n — u n w n . Since |u> n | = 1, we have 

j(v n ) = j(u n ) + p 2 J(w n ) -^U + W = j (3.22) 

in L p (uj) for all p < 2. 

To estimate the energy, we again use the fact that \w n \ = 1 to expand: 

i j a(x)\\7v n \ 2 dx = i / a(z) {|Vw„| 2 + p 2 |Vw„| 2 + j(u„) • j(w n )} dx. 

We claim that the last term is negligible. Indeed, from Step 1, j^jr^ = V$ n , with V<I>„ — > 11/ in C 1 , and 
therefore, 

-^2 / a(x)j{u n ) ■ j{w n ) dx = - / V^ipn ■ p 2 n V$„ da; 

= - j [V^n ' V$„ - (1 - pl)V^ n ■ V$„] cfe 

(l-p 2 )VV„-V$„cfe ^0. 



We calculate, 

lim sup-^G K (v n ;F*) 

= hmsup ^ jf a(x) |i|V U „| 2 + i|Vu,„| 2 - fl, • j(v n ) + |B.| 2 |«„| a + ^(K| 2 - l) 2 + ' ? p 2 n 

< \ jf a(») + ^ j[ a(x) (W + \W\ 2 -B..j + \B*\ 2 + a M^T j dx 
= \( a{x)d\p\ + \ f a(x) (\j-BA 2 + a[X? ^' ? ) dx 



J UJ 



24 
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with = curlB*, where we have used (3.20), (3.5), and (3.22). The completes the proof of the Gamma 
convergence result. 



4. The Obstacle Problem 

We now examine the Gamma limit and its minimizers. Since the limiting functional is two-dimensional, 
we simplify notation somewhat, and denote by x — {x\,X2) <E M 2 , and use j, B 6 R 2 , and F = V x B 
dropping the primes and the asterisks. We also associate to the 1-form j its representation as a vector field 
3 = (ji> J2) € K 2 , and write the Jacobian as a scalar measure, J = x j = \{d Xl j2 — d X2 ji]. 
Using convex duality, we may identify the minimizers of the limiting energy 

h(j\ F) = \ \\a(x) V x + i f a(x)\j - B\ 2 dx 



(obtained in Theorem 1.2) as solutions of a two-obstacle problem for Poisson's equation. The first step is to 
rewrite the minimization problem for I x in terms of a scalar potential function. 

Lemma 4.1. The minimizer of 1^ over Z is attained at j = B — ^V^C; where £ £ Hq(uj) minimizes the 
functional 



Eoo((;F) = 



i \\a(x) (V • i VC + F) + 1 1 ^ I VC| 2 dx, if V • i VC + F € Wl(u), 

+00, otherwise, 



for F = curl B. 

We observe that the Jacobian corresponding to the minimizer is given by: 

Proof. Both Jqo, Eao are convex and lower semicontinuous, so each has a minimizer. Let j G Z. We apply 
the Hodge decomposition from Lemma |3.1| to j — B, to obtain 

j-B = -V^C + V?? + W, 
a 

orthogonal with respect to the inner product (•,•), with C € Hq(ui), r\ G H 1 ^), and W £ H. By the 
orthogonality of the decomposition, 

I oo (j;F) = E oo (0 + ~J a(x){|V?7| 2 + |^| 2 } dx. 

In particular, j minimizes 1^ if and only if the associated £ minimizes -Bqo, and both r\, W = 0. 



Note that for minimizers we conclude that div {a{x){j — B)) = in w, and (j — B) ■ v = on du>. These 
two conditions may also be obtained from the Ginzburg-Landau equations by passing to the thin-film limit 
for minimizers of I E}K - 

Proposition 4.2. Any minimizer £ G Hq(ui) of Erx, is also a minimizer of 

1 



min £oo (it), with E^u) 

|«|<a(*)/2 



a(z) 



(4.2) 
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Proof. We write Foo(C) = *(C) + $(0> with 

i r i 

2 



*(0 



o |VC| 2 , *(C):=2 



art (a 



We then calculate the Legendre transform (conjugate function) of each, with respect to the norm ||u|| 
L il Vu l 2 on H o(u)- Clearly, = = § 1 1 1 1 2 . For $, we have: 



4>*(u) = sup 



sup 
sup 



-Vu- vc- - 

a 2 



u [ -F- V • ^VC ) - 



uv - -\\av\\ m 



a{x) V • -VC + F 



1 

2 

uF. 



a{x) ( V- fl VC + F 



£Crt(a 



uF 



If ||u/a||oo < \-, the bracketed expression above is non-positive for any v, and so the supremum is achieved 
at v = 0. On the other hand, for ||ti/a||oo > \i the bracketed expression is unbounded above. Thus, we 
conclude that 



$*(u) 



' f uF, if II ^11 < i 

Jcj 1 II a Moo — 2 ' 

I + oo , otherwise. 
By the Fenchel-Rockefeller Theorem (see |ET99j or |Bre99] for instance, 

min Erx, = — min (fy*(u) + <&*(—«)) = — min / 



a(x) 



iVul 2 - Fit 



and the minimizers coincide. 



The minimization problem (4.2 1 is the two-obstacle problem, and minimizers u€ K, for 

jueJ#( W ): K*)|<^j 

solve the variational inequality 

/ [Vu • V(u - it) - F (u - u)] > 0, for all vE K. (4.3) 

Minimizers solve the Euler-Lagrange equation away from the coincidence sets, where the pointwise constraint 
is attained, along free boundaries. By Theorem 3.2 in Chapter 1 of |Fri82j . there is a unique, regular minimizer 
to problem ( |4.2| ): 

Proposition 4.3. Let a € C (lj) and F 6 C a (uj), for some a > 0. There exists a unique £ 6 iJ^u;) which 
minimizes Eoo. Moreover, £ € 14 72 ' I? (J7) for all q < oo, and solves 



V- iVC + F = 



on {-2^ <C< ^} 



V-iVC + F>0 on{C>-^} 
V .IVC + F<0 on{C<^} 



(4.4) 
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By the regularity (and the Sobolev embedding) £ g C (w), and thus C satisfies both a Dirichlet and 



Neumann condition on the boundary of the coincidence sets, {ICWI = a(x)}. We also note that by (4.1|, 
the limiting Jacobian for minimizers is thus supported on the two coincidence sets, 

f j- a ( x ) 1 r, { /- a ( x ) 

S- := U e u : C = y'}, S+:=lxeuj: ( = -A_i 

The set S + is thus identified with concentrations of positively oriented vortices (at densities on the order of 
log k) and 5_ supports concentrations of antivortices (at densities on the order of log k) for configurations 
(u,A) with bounded I e>K . 

To illustrate the consequences of Theorems |1.1| and \1.2\ and in particular Proposition |1.3| we revisit two 
examples introduced in our study of the e — > limit in our previous paper [AB GS10] . 

4.1. Vortices and antivortices coexisting 

First we consider a film with uniform thickness and parabolic geometry, f(xi,X2) = x\ + x| — 4, and 
g(xi,X2) — f(xi,X2) + 1, for x' — (xi,X2) € B\ = u) the unit disk. We choose the external field h e 
{ffto&K ) 0, o) , and so for this example, 

~f(x') + g(x>y 



-, ( .r 



a(x') = l, /(*')+ gOO = x f +x l and F = -Hd xi 



-2Hxi. 



Let C,h be the solution of (4.4 1 with parameter H, which we vary according to the external field strength, 
and define 



which solves 





- X\ 


= 


on 




- X\ 


> o 


on 


A£h- 


- Xl 


< o 


on 



4JI Iff/ 

with 6 Hq(Bi). Thus we consider a fixed equation and allow the obstacles to vary as the external field 
strength increases. 

First, observe that if (x\, x%) is a solution, then £,h{xi,X2) '■= —£,h{—xi,X2) is also a solution. By the 
uniqueness of solutions, we deduce that = £h, and hence the solution is odd in x\. 

£,h{xi,X 2 ) = -£h(-XI,X2) 

Furthermore, £h S Hq(uj~) for w~ :=wn {x\ < 0} the left half-disk, and thus solves the two-obstacle 
problem in w~ with Dirichlet boundary condition, 

/ -Vfcr • V(» - > - / ^(fr-u). for all u e 1^, 
X- := | V e iJo 1 ^") : |^)| < i j 

We claim that > in w _ . Indeed, we use v = = max{£#, 0} £ -ftT^ in the above variational inequality. 
With £g = max{— £#,0} we obtain the contradiction 



0< f 1 2 < f Xl Q<0, 

J UJ~ J id 
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Fig. 1. In the first example, vortices first appear for H = 3\/3- For H > 3\/3, positively charged vortices occupy a region in the 
left half-disk, and antivortices concentrate in a symmetrically placed domain in the right half-disk. 



unless = almost everywhere. Thus, £h > in ui~ . Applying the strong maximum principle (see |Fri82| . 
p. 22,) we have £h > in bj~ as claimed. 

By symmetry, we conclude that < in the right half-disk, w + , and thus the coincidence sets (if 
nonempty) S + = S+(H) C cj~ and S- = S-(H) C lu + . By restricting to the half-disks, £h solves a 
one-sided obstacle problem. This problem may be solved explicitly in case the obstacle is not attained: for 
H < 3y/3, the solution for this problem (in polar coordinates) is 

£,h = gKl - f 2 ) cos 6», 

and the coincidence sets S± are empty. When H = 3-^/3, the solution is also given as above, with coincidence 
sets consisting of a single point S± (3V3) = { ( T 4g j 0) } . This value of H marks the lower critical field, the 
smallest value for which vortices appear in global minimizers of the energy. By the comparison principle for 
variational inequalites (see |Fri82j ) . the solutions are monotone decreasing in the half-disk and hence 
the coincidence sets S±{H) form two increasing collections, S±{H) C S±(H) for H < H, symmetrically 
placed across the a^-axis. We illustrate this case in Figure [T] Since the limiting vorticity is given by, 

J = \(-^Ch + Hx x ) - ~ (-A£ H + Xl ) 

we conclude that J < on the set S- and J > on the set S+, indicating the presence of positively oriented 
vortices on S + and antivortices on S_. 



4.2. Vortices accumulating on a ring 

We again choose uj — B\, the unit disk, and take a film of uniform width a{x') — 1 defined by the lower 
surface, 

f(xi,x 2 ) = -x 2 (l - ixj - ^xl 

with g(x±,X2) — f(x%,X2) + 1. Taking an applied field of the form h ex = _ff^|^e*2 gives an effective field 
strength (in polar coordinates), 

F = (4r 2 - 1)H. 
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Making a similar change of variables £h = jj(,H as m the previous example, we obtain a two-obstacle problem 
for £h with a fixed equation but obstacles varying in H, 

A6 f +4r 2 -f = on{-^ < ( H < 577} 
-A£„ + 4r 2 - 1 > on {f H > -^} 
-A£ H + 4r 2 -l<0 on{^< 277} 

The beauty of this problem is that it can be solved explicitly, even after the constraints have been attained, 
and the coincidence sets are determined explicitly as functions of H. Indeed, the general solution of the 
(unconstrained) equation is £ = \r 2 (l — r 2 ) + C\ logr + c 2 , and the constants may be calculated to satisfy 
the obstacle constraints. 
We may then assert: 

(i) When H < 8, there is no coincidence set, the solution does not contact either obstacle. 

(ii) When H = 8, the coincidence set S+ is a circle, with radius r = (The coincidence set S- = 0.) The 
solution is still smooth and is represented in Figure [2] (a). 

(iii) For 8 < H < 16, the coincidence set 5+ expands to an annulus, with the outside radius increasing while 
the inner radius remaining fixed in H: 

with R(H) a strictly increasing function of H. In this interval, the lower obstacle is not attained, so 
S- = 0. This is illustrated in figure [2] (b). 

(iv) At H = 16, the solution contacts the lower obstacle, and SL = {(0,0)}. The approximate value of 
R(W) = 0.8229681606. (See figure |J (c).) 

(v) When H > 16, the coincidence set enlarges on all fronts: 

S + - {p+(H) < v sC R(H)}, 
S- ={r^p-(H)}, 

where < p~{H) < | < p + {H) < 4g. Each function p ± (H),R(H) is monotone, with p ± {H) -» \ and 
R(H) — > 1 as H — > 00. The development of both coincidence sets is illustrated in figure |2j (d) and (e). 
As the Jacobian is supported on the coincidence sets, vortices accumulate both in the annular region 
p + (H) < r < R(H) (with positive orientation), and in the disk r < p~(H) as antivortices (with negative 
flux.) 
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